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Abstract 


This  paper  compares  the  Atlantic  Meridional  Overturning  Circulation  (AMOC)  in  global 
simulations  performed  with  the  depth  coordinate  Parallel  Ocean  Program  (POP)  ocean  model 
and  with  the  HYbrid  Coordinate  Ocean  Model  (HYCOM)  under  different  surface  salinity 
boundary  conditions.  When  forced  by  the  Coordinated  Ocean-ice  Reference  Experiment 
(CORE)  repeat  Normal  Year  Forcing,  HYCOM  develops  internal  inter-annual  variability  during 
the  model  spin-up,  while  POP  does  not  and  HYCOM  is  more  sensitive  to  the  salinity  boundary 
condition  in  the  Southern  Ocean.  Otherwise  the  AMOC  and  related  fields  in  the  two  models  are 
qualitatively  similar,  but  neither  is  able  to  maintain  a  non-trivial  AMOC,  because  of  a  positive 
feedback  that  continually  freshens  the  high-latitude  surface  waters.  However,  with  salinity 
restoring  at  the  ocean  surface  the  AMOC  becomes  progressively  stronger  as  the  piston  velocity  is 
increased.  The  different  restoring  strategies  in  POP  and  HYCOM  cause  differences  in  the 
AMOC  simulation.  The  components  of  the  AMOC  and  closely  related  fields,  including  the 
oceanic  deep  convection,  thermohaline  fluxes,  three  dimensional  currents,  water  mass 
distribution  and  overflows,  are  compared  between  the  models  with  different  salinity  boundary 
conditions.  The  comparison  provides  insights  on  the  models’  response  and  increases  our 
understanding  of  ocean  climate  simulations.  It  also  provides  motivation  for  the  future 
development  of  ocean  climate  models  capable  of  simulating  the  AMOC  more  realistically. 

Keywords:  Atlantic  meridional  overturning  circulation,  global  ocean-ice  modeling,  salinity 
boundary  condition,  isopycnic  coordinate  ocean  model,  depth  coordinate  ocean  model, 
circulation  stability 
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1. 


Introduction 


The  Atlantic  Meridional  Overturning  Circulation  (AMOC)  is  one  of  the  most  important 
components  of  the  global  ocean  circulation.  The  AMOC  here  refers  to  the  circulation  associated 
with  the  oceanic  deep  convection  and  North  Atlantic  Deep  Water  (NADW)  formation  in  the  high 
latitudes  of  the  North  Atlantic.  A  vigorous  AMOC  in  climate  models  transports  a  large  amount 
of  heat  and  salt  northward  in  its  upper  branch,  thereby  maintaining  the  high  surface  temperature 
and  salinity  in  the  northern  North  Atlantic  compared  to  other  oceans  at  the  same  latitudes. 
Failure  to  represent  a  strong  AMOC  would  lead  to  a  simulated  climate  state  very  different  from 
reality,  both  in  terms  of  the  mean  state  and  in  climate  variability  (Manabe  and  Stouffer,  1988; 
Schiller  et  al.,  1997;  Vellinga  and  Wood,  2002;  Stouffer  et  al. ,  2006). 

It  is  challenging  to  obtain  a  realistic  AMOC  in  ocean  climate  numerical  models  because 
the  modeled  circulation  is  highly  sensitive  to  the  thermohaline  forcing  at  the  ocean  surface  and 
its  equilibration  is  extremely  slow  (Rahmstorf,  1997;  Griffies  et  al.,  2009).  A  very  small  bias  in 
the  thermohaline  forcing,  especially  in  the  freshwater  flux,  can  cause  a  persistent  drift  and 
eventual  collapse  of  the  AMOC.  One  factor  is  the  lack  of  direct  feedback  between  sea  surface 
salinity  (SSS)  and  the  freshwater  flux  at  the  ocean  surface  to  control  salinity  drift.  Another  factor 
is  a  strong  positive  feedback  between  the  AMOC  and  northward  salt  transport  in  the  North 
Atlantic  (Stommel,  1961).  A  salinity  anomaly  in  the  high-latitude  North  Atlantic  induced  by  the 
former  mechanism  can  be  amplified  through  the  latter,  leading  to  a  large  salinity  drift  (Paiva  and 
Chassignet,  2001).  In  this  sense,  the  AMOC  is  a  relatively  unstable  system,  in  contrast  to  the 
direct  wind-driven  circulation.  Consequently,  uncoupled  ocean  general  circulation  models 
(OGCMs)  often  require  a  salinity  restoring  to  introduce  a  negative  feedback  to  stabilize  the 
AMOC  (Griffies  et  al.,  2009).  The  simulation  of  the  AMOC  is  more  complex  in  fully  coupled 
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atmosphere-ocean  general  circulation  models  (AOGCMs).  Despite  the  lack  of  salinity  restoring, 
at  least  some  of  the  AOGCMs  that  participated  in  the  Fourth  Assessment  Report  (AR4)  of  IPCC 
were  able  to  maintain  a  stable  vigorous  AMOC  in  their  long-term  control  runs  (Schmittner  et  al., 
2005;  Meehl  et  al.,  2007). 

Most  previous  numerical  investigations  on  the  AMOC  have  been  performed  using  depth 
coordinate  ocean  models  (Bryan,  1986;  Stocker  and  Wright,  1991;  Marotzke  and  Willebrand, 
1991;  Mikolajewicz  and  Maier-Reimer,  1994;  Schmittner  and  Weaver,  2001;  Prange  et  al.,  2003; 
Bryan  et  al.,  2006).  These  originate  from  the  classic  Bryan-Cox  model  developed  in  1960s 
(Bryan,  1969),  and  have  a  long  history  of  application  in  ocean  and  climate  modeling.  Compared 
to  the  depth  coordinate  model,  the  isopycnic  coordinate  model  is  distinctly  different  in  its 
vertical  discretization  of  the  ocean  (Chassignet  et  al.,  1996,  2003;  Griffies  et  al.,  2000;  Bleck, 
2002).  The  depth  coordinate  model  solves  the  equations  using  Eulerian  methodology  in  the 
vertical  direction,  whereas  the  isopycnic  coordinate  model  uses  Lagrangian  (Adcroft  and 
Hallberg,  2006). 

Although  the  development  of  isopycnic  coordinate  models  lagged  depth  coordinate 
models,  this  class  of  models  and  their  derivatives  (i.e.  hybrid  or  quasi-isopycnic  coordinate 
models)  have  matured  to  the  point  that  they  are  routinely  used  in  both  ocean  and  climate 
modeling  (Chassignet  et  al.,  1996,  2003;  Furevik  et  al.,  2003;  Hallberg  and  Gnanadesikan, 
2006;  Sun  and  Bleck,  2006;  Griffies  et  al.,  2009).  Isopycnic  coordinate  ocean  models  preserve 
water  masses  as  they  move  along  isopycnals  and  therefore  more  naturally  represent  the  ocean 
interior.  This  feature  may  be  particularly  advantageous  in  the  simulation  of  the  AMOC 
(Chassignet  et  al.,  2000;  Willebrand  et  al.,  2001).  Considering  that  subgrid-scale 
parameterization  is  a  function  of  the  vertical  coordinate,  it  is  of  great  interest  and  importance  to 
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study  the  performance  of  the  isopycnic  coordinate  model  and  its  derivatives  by  comparing  their 
simulations  of  the  AMOC  with  those  from  the  depth  coordinate  model. 

In  the  present  study,  we  use  the  Parallel  Ocean  Program  (POP)  ocean  model  and  the 
HYbrid  Coordinate  Ocean  Model  (HYCOM)  to  simulate  the  AMOC.  POP  is  a  depth  coordinate 
model  while  HYCOM  is  primarily  an  isopycnic  coordinate  model,  with  the  vertical  coordinate 
smoothly  reverting  to  pressure  coordinate  in  the  mixed  layer  and  unstratified  oceans  (Chassignet 
et  al.,  2006).  Both  models  are  coupled  to  the  Community  Sea  Ice  Model  (CSIM)  and  forced  by 
repeat  Normal  Year  Forcing  (NYF),  as  described  by  Large  and  Yeager  (2004)  and  used  in  the 
Coordinated  Ocean  Research  Experiments  of  Griffies  et  al.  (2009).  In  all  cases  they  are 
initialized  with  the  same  ocean  tracer  climatologies  and  a  state  of  rest.  By  similarly  configuring 
HYCOM  and  POP,  we  investigate  the  impact  of  the  vertical  coordinate  and  salinity  boundary 
conditions  on  the  simulation  of  the  AMOC.  The  goals  are  to  increase  confidence  in  ocean  and 
climate  model  solutions  and  to  advance  ocean  model  development.  The  paper  is  organized  as 
follows:  Section  2  describes  the  models  and  forcing.  Section  3  presents  and  compares  the 
simulation  results  from  HYCOM  and  POP,  followed  by  a  discussion  and  conclusion. 

2.  Description  of  Models  and  Atmospheric  Forcing 

The  depth  coordinate  POP  here  is  version  2,  while  the  one  used  by  Griffies  et  al.  (2009) 
is  version  1  (Smith  and  Gent,  2004)  '.  The  model  employs  Arakawa  B-grid  and  is  configured 
with  the  North  Pole  coordinate  singularity  displaced  over  Greenland.  The  ocean  resolutions  are 
1.125°  in  the  longitudinal  direction,  while  the  latitudinal  resolution  varies  from  0.27°  at  the 
equator  to  0.64°  (totally  320x384  grid  points).  There  are  40  vertical  levels  with  the  thickness  of 

1  POPl.O  and  POP2.0  use  different  advection  schemes.  Some  parameter  values  for  Gent-McWilliams  and 
anisotropic  viscosity  schemes  differ  between  the  two  versions.  They  also  employ  different  salinity  restoring 
strategies  under  sea  ice. 
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the  levels  ranging  from  10  m  at  the  surface  to  250  m  in  the  deep  ocean  (Table  1).  POP  uses  an 
isopycnal  eddy  transport  parameterization  (Gent  and  McWilliams,  1990)  and  an  anistropic 
horizontal  viscosity  (Large  etal.,  2001)  as  generalized  by  Smith  and  McWilliams  (2003).  The  K- 
profile  parameterization  (KPP)  of  vertical  mixing  (Large  et  al.,  1994)  and  an  idealized  diurnal 
cycle  of  solar  forcing  are  also  implemented  in  POP  (Danabasoglu  et  al.,  2006).  The  coefficients 
of  lateral  viscosity  depend  on  the  local  deformation  rate  (Smagorinsky,  1993).  The  background 
vertical  diffusivity  and  viscosity  are  respectively  0. 1  and  1 .0  cm"  s'  .  The  barotropic  equation  is 
solved  using  an  implicit  free-surface  formulation  (Dukowicz  and  Smith,  1994)  and  the 
freshwater  flux  at  the  ocean  surface  is  formulated  as  a  virtual  salt  flux  (Yin  et  al.,  2009). 

HYCOM  (version  2.2)  employs  Arakawa  C-grid  and  is  configured  on  the  same  tracer 
mesh  as  POP.  It  has  32  hybrid  layers  with  the  layer  thickness  ranging  from  3  dbar  (~3  m)  at  the 
surface  to  several  hundred  dbar  in  the  ocean  interior.  Each  coordinate  surface  is  assigned  a 
reference  isopycnal  (Table  2).  In  the  mixed  layer  or  unstratified  oceans,  vertical  grid  points  are 
constrained  to  remain  at  a  fixed  depth,  but  are  allowed  to  join  and  follow  their  reference 
isopycnals  over  the  adjacent  deep  ocean  (Bleck,  2002;  Chassignet  et  al.,  2003,  2006).  Mesoscale 
eddy  stirring  is  implemented  via  layer  interface  smoothing.  Here  the  model  uses  the  standard 
nonslab  KPP  scheme  (Halliwell,  2004).  The  background  vertical  diffusivity  and  viscosity  are  the 
same  as  those  used  by  POP.  As  for  POP,  HYCOM  employs  the  virtual  salt  flux  method.  For 
overflows,  HYCOM  has  a  vertical  resolution  that  naturally  migrates  to  the  density  front  atop  a 
gravity  current  (Hallberg,  2000).  Isopycnic  coordinate  models  have  been  shown  to  produce 
plausible  simulations  of  outflows,  including  the  correct  level  of  entrainment,  although  such 
models  do  show  sensitivity  to  such  key  parameters  as  the  entrainment  rate  and  the  critical 
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Richardson  number  (Xu  et  al.,  2007).  An  improved  representation  of  overflows  in  HYCOM 
could  have  a  significant  impact  on  the  simulations  of  the  AMOC  and  the  mean  climate  state. 

The  prescribed  NYF  surface  forcing  includes  incident  short-wave  insolation  and 
downwelling  long-wave  radiation,  surface  air  temperature  and  humidity,  sea  level  pressure,  wind 
stress,  precipitation  and  runoff  (Large  and  Yeager,  2004).  The  heat  flux  at  the  ocean  surface  is 
determined  by 

Fheat  =  Kad  +  Fsen(SST)  +  Flat(SST)  -  <j(SST)4  ,  (1) 


where  <r  is  the  Stefan-Boltzmann  constant  and  F*ad  is  the  sum  of  the  above  prescribed  radiative 
fluxes  plus  reflected  solar,  with  the  asterisk  indicating  no  year-to-year  variations.  The 
dependence  of  the  upward  long-wave  radiation  and  bulk  calculations  of  both  the  sensible  ( Fsen ) 

and  latent  ( FIat )  heat  fluxes  on  model’s  prognostic  SST  gives  the  flux  inter-annual  variability. 
The  surface  freshwater  water  flux  (not  including  restoring  flux)  is  given  by 

f  =  f*  +f  +  J7* 

water  precip  evap  runoff 

Fevap(SST)  =  Flat/A 


where  X  is  the  latent  heat  of  vaporization  and  both  precipitation  ( F*reci  )  and  runoff  ( F*unoff)  are 
prescribed  with  no  inter-annual  variability.  The  heat  and  water  fluxes  are  partially  coupled 
through  the  calculation  of  evaporation  ( Fevap )  and  latent  heat.  The  freshwater  flux  is  converted  to 

a  virtual  salt  flux  with  the  reference  salinity  ( Sref )  set  to  35  psu: 

F  =  —F  ■ 9  Hi 

salt  water  ^  ref  *  ^  ' 

In  addition,  there  can  be  a  salinity  restoring  term  Frestore  applied  at  the  ocean  surface.  Although 
the  salinity  restoring  is  unphysical,  its  use  in  OGCM  maintains  a  more  realistic  SSS  (Griffies  et 
al.,  2009). 
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The  expression  for  Frestore  differs  between  HYCOM  and  POP.  In  HYCOM,  the  restoring 
salt  flux  is  standard  and  is  calculated  according  to 

FZT  =  ypis,on(SSSohs  -  SSSsJ  ,  (4) 

where  SSSobs  and  SSSsim  are  the  observed  and  simulated  values  respectively,  and  Vpiston  is  a 

piston  velocity.  On  the  other  hand,  POP  uses  anomaly  restoring  which  is  expressed  as, 

pop  _  y  /coo'  _  999'  t 

r restore  *  pistorA^^^  obs  >300  sim) 

=VpislJSSS'obs  - SSSsim(x, y, t)  +  SSS~(t )] 

where  SSS'obs  and  SSS'sjm  are  the  deviations  from  the  global  means  ( SSS )  at  each  time.  The 
anomaly  restoring  in  POP  therefore  ensures  a  zero  net  salt  flux  associated  with  Freslore.  The  salt 

flux  from  anomaly  restoring  depends  on  the  simulated  global  mean  SSS.  It  decreases  with  the 
decrease  in  the  global  mean  SSS,  making  anomaly  restoring  less  effective  than  standard  restoring 
(shown  later).  Therefore,  standard  restoring  is  designed  to  prevent  the  drift  of  the  global  mean 
SSS,  while  allowing  the  drift  of  the  global  mean  salinity.  In  contrast,  anomaly  restoring  is 
designed  to  eliminate  the  effect  of  the  restoring  salt  flux  on  the  global  mean  salinity,  while 
allowing  the  drift  of  the  global  mean  SSS.  Due  to  this  difference  the  restoring  cases  of  HYCOM 
and  POP  are  not  directly  comparable. 

Three  cases  were  integrated  with  HYCOM  and  three  with  POP:  no  salinity  restoring 
( Vpiston  =  0 );  a  weak  salinity  restoring  (Vpiston  =  50m/ 4years);  and  a  strong  salinity  restoring 

( V  iston  =  50m /  30(klays ).  The  restoring  term  is  applied  globally  and  uniformly,  including  under 

sea  ice.  The  comparison  of  the  two  no  restoring  cases  is  used  to  show  that  the  impact  of  the 
different  vertical  coordinate  is  small,  a  conclusion  already  put  forward  by  Griffies  et  al.  (2009). 
The  four  other  restoring  integrations  reveal  a  much  larger  effect  of  the  salinity  boundary 
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condition  on  the  simulation  of  the  AMOC.  HYCOM  is  used  to  look  at  the  impact  of  standard 
restoring  and  POP  to  look  at  the  impact  of  anomaly  restoring.  Due  to  the  computational  expense, 
all  cases  were  run  for  only  150  years,  which  is  sufficient  to  show  differences  in  behavior,  but 
may  not  be  long  enough  for  the  equilibration  of  the  AMOC  in  some  models  (Griffies  et  al., 
2009).  Our  conclusions  in  this  paper  are  based  on  the  150-year  integrations. 

3.  Results 

Figure  1  shows  the  time  evolution  of  the  AMOC  as  simulated  by  HYCOM  and  POP  with 
different  salinity  restoring.  Starting  from  a  resting  state,  the  initial  density  field  almost 
immediately  establishes  the  AMOC.  During  the  first  several  years,  there  is  relatively  little 
sensitivity  of  the  AMOC  to  the  salinity  boundary  condition.  After  its  establishment,  however,  the 
AMOC  cannot  be  sustained  without  a  salinity  restoring.  A  difference  that  may  be  attributable  to 
the  vertical  coordinate  is  the  inter- annual  variability  of  the  AMOC  during  model  spin-up.  With 
the  repeat  annual  cycle  of  NYF,  such  variability  can  only  be  generated  internally  and  is  only  seen 
in  HYCOM.  The  inter-annual  variability  of  the  AMOC  is  quite  pronounced  during  the  first  100 
years  in  the  strong  restoring  case.  The  initial  development  of  the  AMOC  variability  seems 
related  to  the  adjustment  of  the  isopycnic  surfaces  in  HYCOM  (shown  later). 

3.1.  No  Salinity  Restoring  Case 

With  no  salinity  restoring  the  AMOC  rapidly  spins  down  to  less  than  4  Sv  in  both 
HYCOM  and  POP  after  100  years  (Fig.  1).  The  AMOC  index  is  defined  as  the  maximum 
meridional  overturning  streamfunction  value  north  of  30°N  in  the  Atlantic  Ocean.  The 
meridional  streamfunctions  (Figs.  2a  and  2b)  show  the  shutdown  of  the  clockwise  circulation 
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between  about  500  m  and  3000  m,  although  there  is  a  very  weak  remnant  in  POP.  In  HYCOM,  a 
counterclockwise  cell  associated  with  the  circulation  of  the  Antarctic  Intermediate  Water 
(AAIW)  develops  in  the  upper  1000  m  of  the  South  Atlantic  (Fig.  2a),  but  is  confined  to  the 
upper  500  m  in  POP.  It  has  been  shown  that  this  cell  plays  an  important  role  in  causing  multiple 
equilibrium  states  of  the  AMOC  in  some  climate  models  (Yin  and  Stouffer,  2007). 

The  mixed  layer  depths  in  the  high-latitude  North  Atlantic  are  much  too  shallow  in  both 
cases  (Figs.  3b  and  3c)  in  comparison  with  the  observational  estimate  (Fig.  3a),  indicating  a 
complete  cessation  of  the  deep  convection  and  deep  water  formation.  In  both  Figs.  4a  and  4b, 
the  meridional  velocity  in  the  deep  ocean  below  1500  m  is  very  sluggish  without  a  pronounced 
deep  western  boundary  current.  The  Gulf  Stream  is  a  narrow  and  intense  feature,  which  is 
associated  with  strong  titling  and  packing  of  isopycnals.  Its  northward  mass  transport  is 
compensated  by  the  slow,  but  broad  southward  flow  in  the  upper  ocean. 

The  collapse  of  the  AMOC  results  from  the  continuing  decrease  of  the  salinity  in  the 
deep  convection  and  deep  water  formation  region  of  the  North  Atlantic  (Fig.  5b).  In  both 
HYCOM  and  POP,  the  SSS  in  the  region  60°W-30°E,  50°-80°N  decreases  from  about  34.4  psu  to 
a  value  below  31  psu  at  year  150  -  a  roughly  3.5  psu  freshening.  This  freshening  is  sufficiently 
strong  to  shut  down  the  AMOC  in  many  climate  models  (Stouffer  et  al. ,  2006).  In  POP  the 
decrease  of  SSS  is  faster  than  in  HYCOM  during  the  first  70  years,  even  though  the  spin-down 
of  the  AMOC  is  slower  (Fig.  1).  In  terms  of  the  SST  in  the  deep  convection  region,  both  models 
show  a  decrease  from  about  5°C  to  4°C  over  150  years  (Fig.  5a).  Thus,  the  decreases  in  salinity 
and  temperature  contribute  oppositely  to  the  fractional  change  in  the  surface  density  given  by: 

5pt  p0  =  pSS-aST  .  (6) 
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But  with  a  =  1.7xl0-4  “C"1,  /?  =  7.6x10  4  at  these  latitudes,  the  freshening  dominates  the 
cooling  by  an  order  of  magnitude  on  the  right  hand  side  of  Eq.  6,  leading  to  lighter  surface 
water.  The  negative  feedback  between  SST  and  surface  heat  flux  is  evident.  With  the  ocean 
cooling  in  the  deep  convection  region,  the  heat  loss  from  the  ocean  to  the  atmosphere  reduces 
by  roughly  40  W/m  (Fig.  5c),  counteracting  the  decrease  in  SST.  In  contrast,  there  is  no  such 
feedback  between  SSS  and  freshwater  flux,  so  in  both  models  the  freshwater  flux  into  the 
ocean  continually  increases  after  establishment  of  the  AMOC  (Fig.  5d)  in  response  to  the 
reduced  evaporation  due  to  the  cooler  SSTs  (Eqs.  1  and  2),  while  precipitation  remains  as 
prescribed.  The  geographical  distribution  of  the  freshwater  flux  is  quite  similar  between 
HYCOM  and  POP  (Figs.  6a  and  6b). 

With  no  salinity  restoring,  the  geographical  distribution  of  the  SST  biases  compared  to 
the  observation  shows  that  there  is  a  regional  and  moderate  3°C  cooling  in  the  high-latitude 
North  Atlantic  in  both  HYCOM  and  POP  (Figs.  7b  and  7c),  resulting  from  the  shutdown  of  the 
AMOC  and  the  reduction  of  the  northward  heat  transport  in  the  Atlantic.  But  the  spread  of  the 
cooling  to  other  ocean  regions  is  not  evident  due  to  the  constraint  of  the  prescribed  atmosphere. 
In  contrast,  fully  coupled  AOGCMs  showed  that  a  shutdown  of  the  AMOC  can  lead  to  a  10°- 
15°C  ocean  surface  cooling  in  the  northern  North  Atlantic,  and  that  the  cooling  can  spread  over 
the  entire  North  Atlantic  (Yin  and  Stouffer,  2007).  Both  HYCOM  and  POP  show  problems 
common  in  contemporary  coarse  resolution  ocean  climate  models,  such  as  the  displacement  of 
the  Gulf  Stream  and  the  warm  biases  along  the  coastal  region  (Farge  and  Danabasoglu,  2006). 
However,  HYCOM  exhibits  larger  warm  biases  than  in  POP  in  the  Southern  Ocean. 

In  terms  of  SSS,  a  significant  freshening  occurs  in  the  high-latitude  North  Atlantic  in 
both  models  (Figs.  8b  and  8c).  This  large  freshening  illustrates  that  SSS  can  evolve  more  freely 
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than  SST  due  to  the  mixed  boundary  conditions  (Mikolajewicz  and  Maier-Reimer,  1994).  After  a 
freshwater  cap  forms  in  the  northern  North  Atlantic,  it  propagates  southward  with  the  gyre 
circulation  into  the  tropical  and  South  Atlantic.  The  pattern  is  very  similar  to  the  results  from  the 
so-called  “water-hosing”  experiment  (Yin  and  Stouffer,  2007),  although  no  external  freshwater 
addition  is  imposed  in  the  present  study. 

The  most  significant  difference  between  the  formulations  of  HYCOM  and  POP  is  the 
vertical  coordinate.  So  it  is  of  interest  to  compare  the  vertical  distribution  of  temperature,  salinity 
and  therefore  water  masses  as  modeled  by  HYCOM  and  POP.  In  HYCOM,  the  vertical 
coordinate  is  isopycnic  in  the  ocean  interior,  but  smoothly  reverts  to  pressure  coordinate  in  the 
upper  and  polar  oceans,  where  the  ocean  stratification  is  weak  and  deep  convection  takes  place 
(Figs.  9  and  10).  It  should  be  noted  that  due  to  the  shutdown  of  the  AMOC  and  the  surface 
freshening  in  the  high-latitude  North  Atlantic,  the  vertical  coordinate  in  the  upper  ocean  north  of 
40°N  is  isopycnic  in  HYCOM  (Figs.  9b  and  10b).  In  low  and  middle  latitudes,  isotherms  are 
almost  parallel  to  the  isopycnic  interfaces,  as  temperature  dominates  salinity  in  determining  the 
seawater  density  there.  The  upper  warm  water  penetrates  more  deeply  in  POP  (Fig.  9c). 

In  terms  of  salinity,  the  freshening  occurs  only  in  the  surface  layer  in  POP,  so  that  the 
salinity  below  300  m  is  high  when  compared  to  that  in  HYCOM  (Figs.  10b  and  10c).  In  contrast, 
the  freshening  penetrates  more  deeply  in  HYCOM,  so  the  salinity  in  the  entire  North  Atlantic  is 
lower  (Fig.  10b).  In  the  South  Atlantic,  the  northward  intrusion  of  AAIW  is  distinguishable  in 
both  HYCOM  and  POP,  as  shown  by  the  fresh  tongue  in  500-1000  m  depth.  AAIW  reaches 
further  north  in  HYCOM,  so  that  it  connects  to  the  freshening  region  in  the  North  Atlantic. 
However,  the  northward  intrusion  of  AAIW  is  relatively  weak  in  almost  all  model  runs  when 
compared  to  the  observation  (Fig.  10a). 
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3.2.  Salinity  Restoring  Cases 

The  AMOC  and  related  fields  are  much  more  sensitive  to  the  salinity  boundary  condition 
than  to  the  vertical  coordinate’s  choice.  In  all  restoring  cases,  the  upwelling  of  NADW  mainly 
occurs  outside  of  the  Atlantic  basin  and  the  AMOC  is  the  dominant  feature  in  the  Atlantic  (Fig. 
2).  However,  coordinate  differences  may  contribute  to  the  large  Southern  Ocean  changes  in 
ocean  temperature  (Figs.  7  and  9)  and  salinity  (Fig.  10),  seen  in  HYCOM  relative  to  POP. 
Otherwise,  salinity  restoring  can  be  described  as  systematically  increasing  the  strength  of  the 
AMOC  at  150  years,  with  weak  anomaly,  weak  standard,  strong  anomaly  and  strong  standard 
restoring  producing  progressively  stronger  AMOC  (Fig.  2).  This  pattern  of  response  is  also 
indicated  by  the  inter-basin  exchange  across  30°S  and  the  high-latitude  salinities  (Fig.  10). 

Some  responses  are  characterized  by  another  pattern,  from  weak  anomaly  to  strong 
anomaly  to  weak  standard  and  finally  to  strong  standard  restoring.  Notable  examples  of  this 
pattern  are  deep  convection  implied  by  mixed  layer  depths  (Fig.  3),  the  deep  western  boundary 
current  (Fig.  4),  the  SSS  at  years  141-150  (Fig.  5),  and  the  SST  cold  bias  (Fig.  7).  The  overflow 
across  the  Greenland-Iceland-Scotland  ridge  displays  a  more  complicated  pattern.  It  is  evident 
only  in  the  four  restoring  cases  of  Fig.  10  and  characterized  by  the  fresh  tongue  around  50°N  and 
between  1500  to  3000  m.  Further  results  are  now  discussed  separately  for  HYCOM  and  POP, 
since  they  use  different  formulas  for  salinity  restoring  (Eqs.  4  and  5,  respectively). 

3.2.1.  Standard  Restoring  in  HYCOM 

The  AMOC  can  be  sustained  in  HYCOM  once  a  weak  salinity  restoring 
( VptMon =  50m/  4 years)  is  applied  at  the  ocean  surface.  The  AMOC  first  spins  down  with  a 
slower  rate  of  decrease  than  in  the  no  restoring  case.  It  eventually  stabilizes  at  about  12  Sv  after 
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100  years  (Fig.  1).  It  should  be  noted  that  the  definition  of  the  AMOC  index  here  differs  from 
the  one  used  by  Griffies  et  al.  (2009).  Therefore,  Fig.  1  here  cannot  be  compared  directly  to  Fig. 
25  in  Griffies  et  al.  (2009).  Due  to  the  relatively  weak  AMOC,  the  intra-  and  inter-basin 
exchanges  of  seawater  are  also  weak.  There  is  approximately  4  Sv  inter-hemispheric  exchange 
of  seawater  and  8  Sv  of  exchange  at  30°S  associated  with  the  AMOC  (Fig.  2c).  The  AMOC 
remains  active  because  the  fresh  salinity  biases  seen  in  Fig.  8b  are  much  reduced  by  the  salinity 
restoring  (Fig.  8d).  The  SSS  in  the  deep  convection  region  decreases  only  by  0.1  psu  during  the 
150-year  integration,  in  contrast  to  more  than  3  psu  in  the  no  restoring  case  (Fig.  5b).  Thus  the 
salinity  restoring  is  effective  at  modifying  the  total  water  flux  in  the  deep  convection  region 
(Figs.  5d  and  6c)  and  therefore  controlling  salinity  drift  in  HYCOM.  As  a  result  of  the  active 
AMOC,  a  deep  western  boundary  current  is  simulated  just  beneath  the  Gulf  Stream  at  about  2500 
m  depth  (Fig.  4c).  Meanwhile,  the  Gulf  Stream  accelerates,  transporting  more  water  mass 
northward  to  compensate  both  the  upper  southward  gyre  circulation  and  the  deep  western 
boundary  current. 

The  deep  convection  in  the  North  Atlantic  is  realistic  and  takes  place  in  4  regions:  near 
the  south  tip  of  Greenland;  the  Irminger  Sea;  south  of  Iceland;  and  the  Nordic  Seas  (Figs.  3a  and 
3d).  Due  to  the  deep  convection,  upper  ocean  stratification  is  weak  in  the  high-latitude  North 
Atlantic  and  the  vertical  coordinate  in  HYCOM  transitions  from  isopycnals  to  pressure 
coordinate  (Figs.  9d  and  lOd).  The  mode  water  in  the  subtropical  North  Atlantic,  which  is 
characterized  by  salinity  greater  than  35  psu,  can  penetrate  to  a  greater  depth  (Fig.  lOd). 
Generally,  the  spread  of  different  water  masses  along  different  isopycnals  is  more  pronounced  in 
HYCOM  than  in  POP. 
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The  AMOC  is  most  vigorous  and  realistic  in  the  strong  salinity  restoring  case 
( Vpisron =  50/n  /  300days ).  The  mean  strength  of  the  AMOC  during  the  last  20  years  of  the 

integration  is  about  17  Sv  (Fig.  1),  consistent  with  observational  estimates  (Ganachaud,  2003; 
Talley,  2003).  The  intra-  and  inter-basin  exchanges  of  seawater  are  strong,  leading  to  large  heat 
and  freshwater  redistribution  in  the  ocean.  Consistent  with  a  strong  AMOC,  the  speed  of  the 
deep  western  boundary  current  can  exceed  3  cm/s  at  30°N  (Fig.  4e).  A  large  part  of  the 
convection  in  the  Labrador  Sea  and  Nordic  Seas  is  open  ocean  convection  (Fig.  3f).  While  the 
strong  restoring  keeps  the  simulated  SSS  very  close  to  the  observed  (Fig.  8f),  the  SST  biases 
persist  (Figs.  7b,  7d  and  7f),  implying  that  they  are  not  caused  by  the  salinity  drift. 

The  AMOC  in  HYCOM  shows  some  variability  during  the  spin-up  period,  which 
becomes  weaker  after  100  years  (Fig.  1).  Some  fully  coupled  AOGCMs  show  much  stronger 
variability  of  the  AMOC  (Stouffer  el  al.,  2006).  The  mechanism  that  causes  variations  of  the 
AMOC  in  the  uncoupled  OGCM  differs  from  that  in  the  fully  coupled  AOGCM.  The  variability 
of  the  AMOC  in  many  AOGCMs  is  the  ocean  response  to  the  noisy  atmospheric  fluxes 
(Delworth  et  al.,  1993;  Griffies  and  Tziperman,  1995;  Greatbatch  and  Delworth,  2000).  Because 
the  atmosphere  is  prescribed  in  the  present  study,  the  variability  of  the  AMOC  in  HYCOM  is  an 
ocean-only  mode  which  is  induced  by  ocean  internal  processes.  It  is  likely  related  to  the  initial 
adjustment  of  isopycnic  surfaces.  The  thickness  of  the  North  Atlantic  Deep  Water  layers 
(36.38<ci2<36.97,  layers  18-25)  at  30°N  of  the  Atlantic  undergoes  large  variability  during  the 
first  100  years  (Fig.  11).  This  adjustment  of  isopycnic  surfaces  also  occurs  in  other  studies  with 
HYCOM  (Sun  and  Bleck,  2006;  Chassignet  et  al.,  2006).  The  layer  thickness  directly  influences 
the  meridional  mass  transport  in  the  deep  ocean,  contributing  to  the  variability  of  the  AMOC.  In 
addition,  the  mechanism  may  be  linked  to  wave  sloshing  and  needs  to  be  further  investigated. 
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3.2.2.  Anomaly  Restoring  in  POP 

The  AMOC  is  active,  but  weak  (6  Sv)  in  the  weak  salinity  restoring  case  (Fig.  1).  Figure 
2d  shows  that  this  restoring  adds  about  2  Sv  of  overturning  north  of  60°N,  which  gives  most  of 
the  increase  in  the  streamfunction  above  2500  m  in  the  South  Atlantic.  This  picture  is  consistent 
with  NADW  formation  only  at  the  northern  Nordic  Seas,  without  any  deep  convection  and  deep 
water  formation  in  the  Labrador  Sea  and  the  region  south  of  the  Greenland-Scotland  ridge  (Fig. 
3e).  The  POP  results  presented  in  this  paper  are  not  directly  comparable  to  those  presented  in 
Griffies  et  al.  (2009)  because  a)  of  the  use  of  POP  version  2  instead  of  version  1  and  b )  POP  in 
Griffies  et  al.  (2009)  employs  different  salinity  restoring  under  sea  ice.  Weak  anomaly  restoring 
is  not  very  effective  at  modifying  the  water  flux  in  the  deep  convection  region  (Figs.  5d  and  6d) 
and  at  preventing  a  salinity  drift  (Figs.  5b  and  8e).  The  main  reason  is  that  the  global  mean  SSS 
decreases  by  about  0.3  psu  (Fig.  12)  and  the  resulting  restoring  salt  flux  inhibits  the  deep 
convection  in  the  Labrador  Sea  and  south  of  the  Greenland-Scotland  ridge  (Eq.  5  and  Fig.  3e). 
The  SSS  in  the  deep  convection  region  decreases  by  about  1  psu  during  the  first  70  years  before 
the  salinity  restoring  tends  to  prevent  any  further  decrease  (Fig.  5b).  Geographically,  there  is 
still  a  pronounced  freshening  north  of  40°N  by  the  end  of  the  integration  (Figs.  8e  and  lOe).  The 
cooling  of  about  0.5°C  partially  counteracts  the  freshening  and  stabilizes  the  surface  water 
density  (Fig.  5a).  Compared  to  that  in  the  no  restoring  case,  the  reduction  of  the  heat  loss  is 
smaller  due  to  the  smaller  cooling  (Fig.  5c). 

With  strong  anomaly  restoring,  the  AMOC  spins  down  very  slowly  and  reaches  about  14 
Sv  by  the  end  of  the  integration  (Fig.  1).  This  stronger  AMOC  results  from  the  more  realistic 
SSS  in  the  deep  water  formation  region  (Figs.  5b  and  8g).  This  state  is  achieved  while  not 
allowing  restoring  to  change  the  global  ocean  salt  content.  Since  the  global  mean  SSS  decreases 
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by  about  0.2  psu  (Fig.  12),  this  makes  the  strong  anomaly  restoring  less  effective  than  the  strong 
standard  restoring  at  modifying  the  freshwater  flux  (Fig.  6).  The  slow  weakening  of  the  AMOC 
during  the  first  100  years  is  also  a  result  of  a  warming  in  the  deep  convection  region  (Fig.  5a). 
The  heat  loss  to  the  atmosphere  increases  during  the  first  10  years  and  then  gradually  returns  to 
the  initial  value  (Fig.  5c),  leading  to  a  stabilization  of  the  AMOC. 

Vigorous  deep  convection  and  deep  water  formation  takes  place  in  the  Labrador  Sea  and 
Nordic  Seas,  and  in  the  region  south  of  the  Greenland-Scotland  ridge  (Fig.  3g).  Notice  that  the 
deep  convection  sites  in  the  Labrador  Sea  and  Nordic  Seas  are  close  to  the  sea  ice  margin  where 
sea  ice  melting  takes  place.  The  deep  convection  is  mainly  controlled  by  intense  winter  heat  loss. 
In  the  strong  restoring  case,  the  ocean  interior  of  the  northern  North  Atlantic  is  warmer  and  more 
saline  due  to  the  strong  northward  transport  of  heat  and  salt  by  the  AMOC  (Figs.  9  and  10). 

4.  Discussion  and  Conclusion 

The  present  study  uses  two  distinctly  different  ocean  models  to  simulate  the  Atlantic 
meridional  overturning  circulation.  One  model  is  the  depth  coordinate  POP;  the  other  is  the 
quasi-isopycnic  coordinate  HYCOM.  By  comparing  the  simulation  of  the  AMOC  and  closely 
related  fields,  the  present  study  is  the  first  step  to  achieve  our  ultimate  goal,  that  is,  to 
systematically  investigate  the  impact  of  the  vertical  coordinate  and  salinity  boundary  condition 
on  ocean  climate  simulations.  In  general,  the  performances  of  the  two  models  are  qualitatively 
similar  in  terms  of  the  AMOC  modeling.  By  running  similarly  configured  HYCOM  and  POP 
coupled  to  the  same  sea-ice  model  with  identical  prescribed  atmospheric  forcing,  it  is  found  that 
neither  HYCOM  nor  POP  can  maintain  an  active  AMOC  without  the  application  of  salinity 
restoring.  The  two  no  salinity  restoring  cases  are  the  only  two  directly  comparable  cases  in  this 
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study.  It  shows  that  the  two  models  behave  in  a  similar  fashion,  with  all  the  differences  between 
them  being  smaller  than  the  differences  induced  by  salinity  restoring  (Figs.  2-10). 

With  no  salinity  restoring,  the  shutdown  of  the  AMOC  appears  to  be  a  consequence  of 
it  being  established  with  too  little  northward  transport  of  heat  and  salt  into  the  high  latitude 
source  regions  of  the  North  Atlantic  to  balance  the  50  to  60  W/nC  of  surface  heat  loss  and  the  2 
to  3  mm/day  of  freshwater  input.  As  the  SST  and  SSS  decrease  in  response,  there  is  less  net 
heat  loss  and  hence  a  negative  feedback  to  the  SST  decrease,  but  the  reduced  evaporation 
continues  to  add  to  the  freshening.  The  salinity  decrease  and  domination  of  salinity  over 
temperature  in  determining  the  water  density  combine  to  produce  lighter  surface  water,  which 
leads  to  a  weaker  AMOC  and  even  less  northward  transport  of  heat  and  freshwater.  This 
positive  feedback  continues  to  weaken  the  AMOC  until  it  collapses.  After  150  years,  the 
approximate  1°C  decrease  in  SST  reduces  the  heat  loss  by  30  to  40  W/m2.  The  approximate  3 
psu  lower  salinity  means  that  the  density  is  lighter  by  more  than  2  kg/m  .  SSS  and  surface 
density  will  presumably  continue  to  decrease  until  the  gyre  circulation  is  able  to  transport 
sufficient  salt  to  maintain  the  fresher  and  lighter  surface  water. 

Salinity  restoring  is  a  simple,  but  artificial  means  of  overcoming  the  above  positive 
feedback  by  introducing  a  negative  feedback  into  the  system.  The  AMOC  remains  active  once 
salinity  restoring  is  applied,  and  is  strong  in  the  strong  restoring  cases.  Strong  restoring  is  most 
effective  when  applied  in  the  standard  fashion  than  applied  as  anomaly.  The  reason  is  that  the 
anomaly  restoring  flux  depends  on  the  simulated  global  mean  SSS,  which  decreases  in  all 
cases  of  the  POP  simulation.  Griffies  et  al.  (2009)  discussed  other  methods  of  salinity  restoring 
and  present  the  sensitivity  of  other  models.  They  showed  that  the  standard  POP  salinity  boundary 
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condition  of  weak  anomaly  restoring  of  SSS  together  with  a  precipitation  adjustment  produces  an 
acceptable  AMOC  (13  Sv)  in  an  earlier  version  of  POP  (version  1)  driven  by  CORE  NYF. 

Due  to  the  sharply  different  vertical  discretization  between  HYCOM  and  POP,  the 
similarities  between  HYCOM  and  POP  indicate  robust  results.  The  results  confirm  that  the 
thermohaline  forcing,  especially  the  salinity  boundary  condition  and  freshwater  forcing  at  the 
ocean  surface,  plays  a  distinctive  and  crucial  role  in  simulating  the  AMOC.  Except  for  the 
salinity  boundary  condition,  all  other  forcings,  including  wind  stress  and  model  parameters  such 
as  for  diapycnal  mixing,  are  the  same  between  different  cases.  The  simulated  anemic,  weak  and 
strong  AMOCs  in  the  three  cases  suggest  that  proper  thermohaline  fluxes  are  among  the  most 
important  and  difficult  issues  in  the  AMOC  modeling,  especially  when  the  ocean  model  is 
coupled  to  the  atmospheric  model  and  when  salinity  restoring  should  not  be  used. 

With  repeat  NYF,  the  AMOC  collapses  in  the  uncoupled  ocean  model  unless  there  is 
some  unphysical  freshwater  flux.  However,  the  AMOC  is  typically  strong  when  ocean  models 
such  as  POP  are  coupled  to  the  atmospheric  models,  even  though  the  atmospheric  fluxes  in  the 
coupled  system  may  deviate  greatly  from  observational  estimates  and  there  is  no  salinity 
restoring.  One  possible  explanation  is  that  the  maintenance  of  the  AMOC  may  depend  critically 
on  extremely  cold  winters,  which  are  missing  in  repeat  NYF.  A  better  understanding  of  the 
reason  for  AMOC  maintenance  in  AOGCMs  would  be  very  useful  for  model  development  and 
ocean  climate  simulation,  where  an  important  goal  must  be  generate  a  stable  and  vigorous 
AMOC  without  any  unphysical  salinity  restoring.  Until  then  it  is  difficult  to  go  beyond  the 
present  study  where  the  simulations  of  the  components  of  the  AMOC,  including  ocean  currents, 
overflow,  and  water  mass  distribution  are  either  far  from  nature,  or  confused  by  unphysical 
boundary  conditions. 
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Figure  Legends 

Figure  1.  Time  evolution  of  the  AMOC  in  the  runs  without  salinity  restoring  (“NO”)  and  with 
weak  (“WEAK”)  and  strong  (“STRONG”)  salinity  restoring.  The  AMOC  index  is  defined  as  the 
maximum  meridional  overturning  streamfunction  value  north  of  30°N  in  the  Atlantic. 

Figure  2.  The  AMOC  streamfunction  patterns  (Sv)  simulated  by  HYCOM  and  POP  with 
different  salinity  boundary  conditions  (mean  of  the  last  10  years).  Positive  values  indicate 
clockwise  circulation,  (a),  (c)  and  (e):  HYCOM;  (b),  (d)  and  (f):  POP;  (a)  and  (b):  no  salinity 
restoring  case;  (c)  and  (d):  weak  salinity  restoring  case;  (e)  and  (f):  strong  salinity  restoring  case. 

Figure  3.  Annual  mean  mixed  layer  depth  (m)  in  the  deep  convection  region  of  the  AMOC.  (a) 
The  observational  estimate;  (b),  (d)  and  (f):  HYCOM;  (c),  (e)  and  (g):  POP;  (b)  and  (c):  no 
salinity  restoring  case;  (d)  and  (e):  weak  salinity  restoring  case;  (f)  and  (g):  strong  salinity 
restoring  case.  The  observational  estimate  of  mixed  layer  depth  is  derived  from  monthly  mean 
climatology  of  potential  temperature  and  salinity  (Griffies  et  al.,  2009). 

Figure  4.  The  meridional  velocity  (cm/s)  at  30°N  of  the  Atlantic  (mean  of  the  last  10  years).  The 
data  are  plotted  on  the  native  grids  of  HYCOM  and  POP.  The  interfaces  of  the  hybrid  layers  in 
HYCOM  are  indicated  by  the  black  curves.  The  numbers  show  the  interface  index. 

Figure  5.  Time  evolution  of  the  SST,  SSS,  heat  and  water  fluxes  averaged  in  the  deep 
convection  and  deep  water  formation  region  of  the  AMOC  (50°-80°N,  60°W-30°E).  Positive 
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fluxes  indicate  a  transport  from  the  atmosphere  to  the  ocean,  (a)  SST  (°C);  (b)  SSS  (psu);  (c)  heat 
flux  (W/m  );  (d)  water  flux  (mm/day). 

Figure  6.  Geographical  distribution  of  the  total  freshwater  flux  (mm/day)  including  the  restoring 
part  in  HYCOM  and  POP. 

Figure  7.  Geographical  distribution  of  the  SST  biases  (°C).  (a)  The  SST  observation  (Steele  et 
al.,  2001);  (b),  (d)  and  (f):  biases  in  HYCOM;  (c),  (e)  and  (g):  biases  in  POP;  (b)  and  (c):  no 
salinity  restoring  case;  (d)  and  (e):  weak  salinity  restoring  case;  (f)  and  (g):  strong  salinity 
restoring  case.  The  values  show  the  simulation  (mean  of  the  last  10  years)  minus  the  observation. 

Figure  8.  Same  as  Fig.  7,  but  for  the  SSS  biases  (psu ).  Notice  that  (d)  and  (e)  use  different 
scales. 

Figure  9.  The  latitude-depth  cross-section  of  ocean  temperature  (°C)  along  30°W  in  the  Atlantic. 
The  data  (shading)  are  plotted  on  the  native  grid  of  HYCOM  and  POP.  Green  lines  show  the 
temperature  biases  (interval=2.0)  and  black  ones  show  sample  isopycnals. 

Figure  10.  Same  as  Fig.  9,  but  for  salinity  (psu).  Green  lines  show  the  salinity  biases 
(interval=0.5). 
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Figure  11.  Changes  in  layer  thickness  (m)  at  30°N  of  the  Atlantic  in  HYCOM.  The  values 
indicate  the  deviation  from  the  long-tenn  mean.  Red  curves  show  the  layers  that  are  occupied  by 
North  Atlantic  Deep  Water. 

Figure  12.  Time  series  of  the  global  mean  SSS  (psu )  in  HYCOM  and  POP. 
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Figure  6 
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Figure  7 
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Figure  8 


Figure  9 
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Figure  10 
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Figure  12 


